
  ;;......................................... 
  ;; Input data 
  ;;......................................... 
  fna = "sgp_2010_cam5_v2.nc"
  fnb = "sgp_20100401_shaocheng.nc"

  fla = addfile(fna,"w")
  flb = addfile(fnb,"r")

  ;;......................................... 
  ;; Base date
  ;;......................................... 
  bdate = flb->bdate
  bdate = (/ 20091201 /)
  fla->bdate=bdate

  ;;......................................... 
  ;; Lat & Lon
  ;;......................................... 
  lat = 36.605
  lat!0 = "lat"
  lat&lat = 36.605

  lon = -97.485
  lon!0 = "lon"
  lon&lon = -97.485

  fla->lon=lon
  fla->lat=lat
  time = fla->time

  lev = flb->lev
  fla->lev=lev

  nt = dimsizes(time)

  tsec = new(nt,"double")
  tsec!0 = "time"
  tsec&time = time
  tsec@long_name = "time" ;
  tsec@units = "seconds" ;

  tsec = (/ fspan(0.,3600.*nt,nt)   /)

  fla->tsec=tsec

